Stochastic background of relic scalar gravitational waves from scalar-tensor gravity 
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A stochastic background of relic gravitational waves is achieved by the so called adiabatically- 
amplified zero-point fluctuations process derived from early inflation. In principle, it provides a 
distinctive spectrum of relic gravitational waves. In the framework of scalar-tensor gravity, we discuss 
the scalar modes of gravitational waves and the primordial production of this scalar component which 
is generated beside tensorial one. We discuss also the upper limit for such a relic scalar component 
with respect to the WMAP constraints. 
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The design and the construction of a number of sensitive detectors for gravitational waves (GWs) is underway today. 
There are some laser interferometers like the VIRGO detector, built in Cascina, near Pisa, Italy, by a joint Italian- 
French collaboration, the GEO 600 detector built in Hanover, Germany, by a joint Anglo-German collaboration, the 
two LIGO detectors built in the United States (one in Hanford, Washington and the other in Livingston, Louisiana) 
by a joint Caltech-Mit collaboration, and the TAMA 300 detector, in Tokyo, Japan. Many bar detectors are currently 
in operation too, and several interferometers and bars are in a phase of planning and proposal stages (for the current 
status of gravitational waves experiments see [H, Q)- The results of these detectors will have a fundamental impact on 
astrophysics and gravitation physics. There will be a huge amount of experimental data to be analyzed, and theorists 
would face new aspects of physics from such a data stream. Furthermore, GW detectors will be of fundamental 
importance to probe the General Relativity or every alternative theory of gravitation [f| • A possible target of 

these experiments is the so called stochastic background of gravitational waves @> fl SllOjll]] • The production of the 
primordial part of this stochastic background (relic GWs) is well known in literature starting from the works by Q 
and which, using the so called adiabatically-amplified zero-point fluctuations process, have shown in two different 
ways how the inflationary scenario for the early universe can, in principle, provide the signature for the spectrum of 
relic GWs. 

In this letter, after a short summary of scalar-tensor gravity admitting the existence of GW scalar modes, the 
primordial production of this scalar component is discussed. Besides, the results are confronted to the WMAP data 
|12l. [l3| in order to achieve an upper limit for the scalar part of relic GWs. 

Scalar-tensor gravity theories are a particular case of Extended Theories of Gravity which are revealing a useful 
paradigm to deal with several problems in cosmology, astrophysics and fundamental physics (for a comprehensive 
discussion see, for example [3, EB|)- In the most general case, considering only the Ricci scalar among the curvature 
invariants, they arises from the action 
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where F is an unspecified function of curvature invariants and of a scalar field and □ is the D'Alembert operator. 
The term C m is the minimally coupled ordinary matter contribution. Scalar-tensor gravity, is recovered from (U) 
through the choice 



F(R,ct>)=f( < p)R-V( ( f>), 
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Considering ((2|), a general action for scalar-tensor gravity in four dimensions is 



S = d A x^J —g 



(3) 



2 



which can be recast in a Brans-Dicke-like form by 

tp = f(4>) , w{<p) = 4^ , w(<p) = vMtp)) 



and then 
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S = J d X\J—g 

By varying the action ([5]) with respect to g M „, we obtain the field equations 

while the variation with respect to tp gives the Klein - Gordon equation 
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We are assuming physical units G = 1, c = 1 and h = 1. Tu™^ is the matter stress-energy tensor and G is a 



dimensional, strictly positive, gravitational coupling constant [1, [||. The Newton constant is replaced by the effective 
coupling 
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eff 
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which is, in general, different from G. General Relativity is recovered for 

1 
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In order to study gravitational waves, we assume first-order, small perturbations in vacuum (T#° = 0). This means 

SV = + h^v , p = p Q + Sp (10) 

and 
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for the self-interacting, scalar-field potential. These assumptions allow to derive the "linearized" curvature invariants 
Rjxvpu , Rfiv and R which correspond to R^ vpa , R^v and R, and then the linearized field equations @, [13] 
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The case u> = const and W = has been analyzed i n [p| considering the so-called "canonical" linearization [IT]]. In 
particular, the transverse-traceless (TT) gauge (see (l7fl ) can be generalized to scalar-tensor gravity obtaining the 
total perturbation of a GW incoming in the z+ direction in this gauge as 

V(* - *) = A + {t - z)e<+> + A x (t - «) e M + $(* - z) e W (14) 

The term A + (t — z)e^ + A x (t — z)e$ describes the two standard (i.e. tensorial) polarizations of a gravitational 
wave arising from General Relativity in the TT gauge [T3|, while the term $(i — z)e p s } is the extension of the TT 
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gauge to the scalar case. This means that, in scalar-tensor gravity, the scalar field generates a third component for 
the tensor polarization of GWs. This is because three different degrees of freedom are present (see Eq.(32) of 0), 
while only two are present in standard General Relativity. 

Let us now take into account the primordial physical process which gave rise to a characteristic spectrum fl sgw for 
the early stochastic background of relic scalar GWs. The production physical process has been analyzed, for example, 
in 0, 0, [HI but only for the first two tensorial components of eq. lfT4|) due to standard General Relativity. Actually 
the process can be improved considering also the third scalar-tensor component. 

Before starting with the analysis, it has to be emphasized that, considering a stochastic background of scalar GWs, 
it can be described in terms of the scalar field $ and characterized by a dimensionless spectrum (see the analogous 
definition for tensorial waves in 0, SUB]) 
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is the (actual) critical energy density of the Universe, H the today observed Hubble expansion rate, and dp sgw is the 
energy density of the scalar part of the gravitational radiation contained in the frequency range / to / + df. We are 
considering now standard units. 

The existence of a relic stochastic background of scalar GWs is a consequence of general assumptions. Essentially it 
derives from basic principles of Quantum Field Theory and General Relativity. The strong variations of gravitational 
field in the early Universe amplifies the zero-point quantum fluctuations and produces relic GWs. It is well known 
that the detection of relic GWs is the only way to learn about the evolution of the very early universe, up to the 
bounds of the Planck epoch and the initial singularity 0, d, 0, [T3, EH ■ It is very important to stress the unavoidable 
and fundamental character of such a mechanism. It directly derives from the inflationary scenario [l3,[23|, which well 
fit the WMAP data in particular good agreement with almost exponential inflation and spectral index « 1, [H, 

A remarkable fact about the inflationary scenario is that it contains a natural mechanism which gives rise to 
perturbations for any field. It is important for our aims that such a mechanism provides also a distinctive spectrum 
for relic scalar GWs. These perturbations in inflationary cosmology arise from the most basic quantum mechanical 
effect: the uncertainty principle. In this way, the spectrum of relic GWs that we could detect today is nothing else 
but the adiabatically- amplified zero-point fluctuations 0, The calculation for a simple inflationary model can be 
performed for the scalar field component of eq. lfl4|) . Let us assume that the early Universe is described an inflationary 
de Sitter phase emerging in a radiation dominated phase 0, [1, [13] • The conformal metric element is 



ds 2 = a 2 (f])[-dT] 2 + dl? 2 + h^(r ] ,^)dxi J -dx v }, 
where, for a purely scalar GW the metric perturbation (fLi)) reduces to 

Following 0, [13], in the de Sitter phase, we have: 
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while, in the radiation dominated phase we have, respectively, 
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t]i is the inflation-radiation transition conformal time and t] Q is the value of conformal time today. If we express the 
scale factor in terms of comoving time cdt = a(t)drj, we have 



a(t) ex enp(Hd s t) , a(t) oc Vt (19) 

for the de Sitter and radiation phases respectively. In order to solve the horizon and flatness problems, the condition 
— 7"^r > 10 27 has to be satisfied. The relic scalar-tensor GWs are the weak perturbations huvij], it) of the metric lfl~8l) 
which can be written in the form 

hfj,u = e$(k)X(rj) exp(t ■ it), (20) 
in terms of the conformal time r] where A: is a constant wavevector. From Eq. ((20|l . the scalar component is 

$(r),le, it) = X(rj)cxp(t ■ it). (21) 
Assuming Y(rf) — a(rj)X(r]), from the Klein-Gordon equation in the FRW metric, one gets 

Y" + (\t\ 2 -— )Y = (22) 
a 

where the prime ' denotes derivative with respect to the conformal time. The solutions of Eq. |22|) can be expressed 
in terms of Hankel functions in both the inflationary and radiation dominated eras, that is: 
For 77 < 771 

X(v) = ^[l + fTd^lexp-ifcfa-Tfc), (23) 

for r\ > i]i 

X{v) = t~v \ a ex P —ik(r] — m) + /3exp ik(n — m), (24) 
a{r)) 

where to — ck/a is the angular frequency of the wave (which is function of the time being k = | k \ constant), a and 
are time-independent constants which we can obtain demanding that both X and dX/dr] are continuous at the 
boundary 77 = 771 between the inflationary and the radiation dominated eras. By this constraint, we obtain 



.y/HdsHa H ds H a H ds H 
a = 1 + l — 2^' /3= ^" (25) 

In Eqs. |25ll . uo — ck/a(r]o) is the angular frequency as observed today, H = c/rj is the Hubble expansion rate as 
observed today. Such calculations are referred in literature as the Bogoliubov coefficient methods 0, E3|- 

In an inflationary scenario, every classical or macroscopic perturbation is damped out by the inflation, i.e. the 
minimum allowed level of fluctuations is that required by the uncertainty principle. The solution |23|) corresponds to a 
de Sitter vacuum state. If the period of inflation is long enough, the today observable properties of the Universe should 
be indistinguishable from the properties of a Universe started in the de Sitter vacuum state. In the radiation dominated 
phase, the eigenmodes which describe particles are the coefficients of a, and these which describe antiparticles are 
the coefficients of /3 (see also [ll|). Thus, the number of particles created at angular frequency u) in the radiation 
dominated phase is given by 

^H/3.| 2 =(^) 2 . (26) 

Now it is possible to write an expression for the energy density of the stochastic scalar relic gravitons background in 
the frequency interval (u>, uj + doj) as 

do 2^ (^^N hH l H o d f (27] 

d Psgw - 2hu [^&)N U - - 4^-3-— - ~l[^rj > (27) 
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Figure 1: The spectrum of relic scalar GWs in inflationary models is flat over a wide range of frequencies. The horizontal axis 
is log 10 of frequency, in Hz. The vertical axis is log 10 £l gs w The inflationary spectrum rises quickly at low frequencies (wave 
which re-entered in the Hubble sphere after the Universe became matter dominated) and falls off above the (appropriately 
redshifted) frequency scale f ma x associated with the fastest characteristic time of the phase transition at the end of inflation. 
The amplitude of the flat region depends only on the energy density during the inflationary stage; we have chosen the largest 
amplitude consistent with the WMAP constrains on scalar perturbations. This means that at LIGO and LISA frequencies, 



< 2.3 * 10' 



where /, as above, is the frequency in standard comoving time. Eq. (|27|) can be rewritten in terms of the today and 
de Sitter value of energy density being 

8nG Pc _ 8TrG Pds 

Ho = ^c^> Hds '^^- (28) 

c 7 

Introducing the Planck density ppianck = tt^ the spectrum is given by 

m _ 1 dp sgw _ f dp sgw _ 16 p ds 

p c d 111 / p c dj 9 ppianck 

At this point, some comments are in order. First of all, such a calculation works for a simplified model that does not 
include the matter dominated era. If also such an era is also included, the redshift at equivalence epoch has to be 
considered. Taking into account also results in we get 

fW/) - ^-^-(1 + Ze q )~\ (30) 
y PPlanck 

for the waves which, at the epoch in which the Universe becomes matter dominated, have a frequency higher than 
H eq , the Hubble parameter at equivalence. This situation corresponds to frequencies / > (1 + Zeq) 1 / 2 !!^. The redshift 
correction in Eq. (f30|l is needed since the today observed Hubble parameter Hq would result different without a matter 
dominated contribution. At lower frequencies, the spectrum is given by UM 

CW/)oc.T 2 . (31) 



As a further consideration, let us note that the results |29|l and |30j) . which are not frequency dependent, does not 
work correctly in all the range of physical frequencies. For waves with frequencies less than today observed Hq, the 
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notion of energy density has no sense, since the wavelength becomes longer than the Hubble scale of the Universe. 
In analogous way, at high frequencies, there is a maximal frequency above which the spectrum rapidly drops to 
zero. In the above calculation, the simple assumption that the phase transition from the inflationary to the radiation 
dominated epoch is instantaneous has been made. In the physical Universe, this process occurs over some time scale 
At, being 

_ a(ti) 1 

J max — J7~ \~\ ) (32 J 

a(t ) At 

which is the redshifted rate of the transition. In any case, £l sg w drops rapidly. The two cutoffs at low and high 
frequencies for the spectrum guarantee that the total energy density of the relic scalar gravitons is finite. For GUT 
energy-scale inflation it is of the order 

Pds « 1CT 12 . (33) 



PPlanc 

These results can be quantitatively constrained considering the recent WMAP release. In fact, it is well known that 
WMAP observations put strongly severe restrictions on the spectrum. In Fig. Q] the spectrum fl sgw is mapped : 
considering the ratio Pds/ ppianck, the relic scalar GW spectrum seems consistent with the WMAP constraints on 
scalar perturbations. Nevertheless, since the spectrum falls off oc f~ 2 at low frequencies, this means that today, at 
LIGO-VIRGO and LISA frequencies (indicated in fig. [1]), one gets 

n sg w(f)h 2 W0 < 2.3 x 1CT 12 . (34) 

It is interesting to calculate the corresponding strain at ~ 100Hz, where interferometers like VIRGO and LIGO reach 
a maximum in sensitivity. The well known equation for the characteristic amplitude @, [l(J [ill adapted to the scalar 
component of GWs can be used: 



1Hz. 



<M/) ~ 1.26 x lO-^(--)^hl Qo n sgw (f), (35) 
and then we obtain 

$ c (100Hz) < 2 x 1CT 26 . (36) 

Then, since we expect a sensitivity of the order of 10 -22 for the above interferometers at w 100Hz, we need to 
gain four order of magnitude. Let us analyze the situation also at smaller frequencies. The sensitivity of the VIRGO 
interferometer is of the order of 10 -21 at w 10 H z and in that case it is 

§ c {100Hz) < 2 x 1CT 25 . (37) 

The sensitivity of the LISA interferometer will be of the order of 10 -22 at w 10~ 3 Hz and in that case it is 

$ c (100Hz) < 2 x 1CT 21 . (38) 

This means that a stochastic background of relic scalar GWs could be, in principle, detected by the LISA interferom- 
eter. 

In summary, the above results point out that a further scalar component of GWs, coming from Extended Theories 
of Gravity, should be seriously considered in the signal detection of interferometers. As already discussed in (_!]], this 
fact could constitute an independent test for alternative theories of gravity or a further probe for General Relativity 
capable of ruling out other theories. 
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